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We address the problem of evaluating the number Sn (t) of distinct sites visited up to time t by N 
noninteracting random walkers all initially placed on one site of a deterministic fractal lattice. For a 
wide class of fractals, of which the Sierpinski gasket is a typical example, wejaropose that, after the 
short-time compact regime and for large N, SN(t) ~ S'jv(t)(l — A), where 5jv(f) is the number of 
sites inside a hypersphere of radius R[ln(N) jc\ , R is the root-mean-square displacement of a single 
random walker, and u and c determine how fast 1 — r t (r) (the probability that site r has been visited 
by a single random walker by time t) decays for large values of r/R: 1 — T t (r) ~ exp[— c(r/R) u ]. 
For the deterministic fractals considered in this paper, u = d m /(d w — 1), d w being the random walk 
dimension. The corrective term A is expressed as a series in ln~ n (iV) ln m ln(iV) (with n > 1 and 
< m < n), which is given explicitly up to n — 2. Numerical simulations on the Sierpinski gasket 
show reasonable agreement with the analytical expressions. The corrective term A contributes 
substantially to the final value of Sjv(i) even for relatively large values of N. 

PACS numbers: 05.40.Fb,05.60.Cd,66.30.Dn 



I. INTRODUCTION 



Random walk theory is a branch of statistical physics with many applications pL|^| . Problems related to a single 
random walker have traditionally been the subject of thorough study, but their generalizations to N > 1 random 
walkers have attracted much less attention, although there are some, generally very recent, exceptions These 
multiparticlc diffusion problems are characterized by the impossibility of being analyzed in terms of the single random 
walker theory, i.e., they can not be solved through simple averaging over the properties of a single random walker, 
even in the noninteracting case. The recent development of experimental techniques allowing the observation of 
events caused by single particles of an ensemble should give additional encouragement to the study of these of 
multiparticle diffusion problems. 

The problem that is the subject of this paper, namely, the evaluation of the average number S]\r(t) of distinct 
sites visited (or territory covered) by N random walkers up to time t, all moving from the same starting site, is a 
clear example of a diffusion problem that can not be solved, or even approximated, from the solution for N = 1, 
Si(t). Even for independent random walkers, the overlap of the regions explored by different walkers prohibits a 
decomposition of Sjsr(t) into single-particle contributions. The origin of the problem of evaluating Sjv(i) is usually 
traced back to the Proceedings of the Second Berkeley Symposium on Mathematical Statistics and Probability where 
the case N = 1 was posed by Dvoretzky and Erdos (TlJ. Since then, the quantity Si(t) has been studied in detail 
and is discussed in general references [l||T^]. For fractal substrates (in particular, for the two-dimensional Sierpinski 
gasket and the two-dimensional percolation cluster at criticality) this problem was first studied by Rammal et al. 
P,[l3|. More recently, Larralde et al. |fl and Havlin et al. |1 studied the problem of evaluating S^(t) when N ^> 1 
noninteracting random walkers diffuse in Euclidean and fractal media, respectively. For fractal lattices with spectral 
dimension d s = 2d f /d w < 2, it was argued § that S N (t) ~ t dl for t < t < t x - In N and S N {t) - t d ^ 2 ln d//tl (N) for 
t x <C t, where d[ = df /d m ; n is the chemical dimension (or topological distance dimension), df is the fractal dimension, 
dmin is the fractal dimension of the shortest path on the fractal, u — d w /(d w — 1) and d w is the diffusion exponent (or 
fractal dimension of the random walk) However, Drager and Klafter |l using scaling arguments have recently 

proposed that 

S N (t) ~ t d °! 2 (\nN) d i! v (1) 

for t x <C t, where v = <4,/(<4; — 1) an d d l w = d w /d min is the chemical-diffusion exponent p|,^i|. Of course, the 
two predictions agree for the media considered, such as Sierpinski gaskets, in the present paper for which d m ; n = 1. 
Therefore, in what follows, we will write dj , u and d w instead of de, v and d l w , respectively. Fractals with <i m i n ^ 1 
will be discussed in the following paper |lq|. 
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As stated above, two time regimes are observed in Sn^): an extremely short-time regime or regime I and a long- 
time regime or regime II separated by the crossover time t x ~ IniV. A further long-time regime, or regime III, is 
observed in Euclidean lattices when the movement of the independent walkers are very far from each other so that 
their trails (almost) never overlap and Sjv'(t) ~ NSi(i) |^,@J|- In the one-dimensional lattice and fractal lattices 
with d s < 2, the trails of the random walkers partially overlap at all times and regime III is never reached. Such is 
the case in this paper where we are concerned only with fractals in which d s < 2. The transition from regime I to 
regime II is easy to understand. In regime I, we have so many particles at every site that all nearest neighbors of the 
already visited sites are necessarily reached by some walker at the next time step. The minimum path length between 
two sites on a fractal is the chemical distance, £, and it is clear that after a time t the visited zone is an hypersphere 
of chemical radius £ = t and, consequently, Sjf (t) ~ t de according to the definition of the chemical exponent [p|,p^[. 
This behavior lasts until the average number of particles per site is of the order of unity, then overlapping is only 
partial and a new regime is established. If z is the coordination number of the lattice (or the average coordination 
number in the case of stochastic fractals), the number of random walkers at every site decreases as N/z l in the very 
short-time regime because the particles are distributed between z sites at every step. The breaking of the overlapping 
regime thus takes place when N/z l ~ 1 or, equivalently, t x ~ InN. 

Regime II is far more interesting and difficult to analyze than regime I due to the nontrivial interplay of the walkers 
in their exploration of the lattice, which leads to a more complex dependence of Sisr(t) on time t and number of 
particles N. In some recent work we have shown that for independent random walks on Euclidean lattices 

there exist important asymptotic corrections to the main term that can not be ignored even for very large number 
of particles as these corrections decay only logarithmically. We will see that this also holds for the two-dimensional 
Sierpinski gasket. 

It should be noticed that, except for the Sierpinski gasket in two dimensions when N — 1 B3|, there has never 
been any discussion about Sjsr(t) focused on deterministic fractals, whether theoretically or numerically. Certainly, 
a dependence on t and N of the main asymptotic term of Siv(t) for large N has been proposed [see Eq. (jl])], but 
nothing is known about the value of its amplitude or prefactor and on the relevance (if any) of the other (corrective) 
asymptotic terms. In this paper we present a procedure for obtaining, for a certain class of fractals, the complete 
asymptotic series expansion of SW(i) when N 3> 1. The procedure gives the main asymptotic term in full, and 
determines the functional form of the corrective terms, which we calculate explicitly up to second order. The fractals 
that we consider in this paper have to satisfy two conditions. First, the number of sites (or volume) V(r) of the fractal 
inside a hypersphere of radius r should be given by 



V(r) = V r d f (2) 

where Vb is a constant characteristic of the fractal substrate; and, second, the probability T t (r) that a site r has not 
been visited by a single random walker by time t should decay for £ = |r|/i? = r/R^S> 1 as 

T t (r) a 1 - AC' m exp(-cO (l + + ■■■), (3) 

where R 2 = 2Dt 2 / d ™ is the mean-square displacement of a single random walker and D is the diffusion constant. It 
should be clear at this point that the above two conditions can only be satisfied approximately: first, because Vq is not 
strictly constant (it exhibits log-periodic oscillations of small amplitude, see Sec. |l|); second, because T t (r) does not 
solely depend on the distance r but also (in general) on the actual location r on the lattice; and, third, because T t (r) 
is not continuos (this fact can be clearly seen in figure 1 of Ref. jl6|]) so that Eq. (||) can only be an approximation 
to the true distribution. The fluctuations in 5jv(i) associated to these effects arc thus not included in our theoretical 
discussion. However, their importance can be gauged by resorting to simulation. For the two-dimensional Sierpinski 
gasket, we found that these fluctuations are indeed relevant and that they can be explained to a large extent as a 
consequence of the log-periodic oscillations of Vq. Finally, there is another difficulty regarding the value of T t (r): 
while its dominant term exp(— c£ u ) is reasonably well established, the value (and even the form) of its subdominant 
factors A^~^ u , etc. is unknown (although an educated guess can be made; see Sec. 0). This means that we 

can be fairly sure of the value of the main term of Si\r(t) because, as we will show, it depends only on the dominant 
term of r t (r). However, the true value of the corrective terms of Sjy(i) is more uncertain as they also depend on 
the subdominant factors of r t (r). Nevertheless, we will see that reasonable choices of values for these subdominant 
factors lead to significant improvements in the estimate of 5/v(i). 

The simulation results for the two-dimensional Sierpinski gasket will reveal the great importance of the corrective 
terms (for example, they account for some thirty per cent of the total value of SW(t) even for N as large as 10 6 ). This 
is not strange because, as we will show in Sec. O, they decay only logarithmically with N. An important consequence 



that we will address in the following paper 
based on "collapsing" the numerical data J 



is that the corrective terms must be taken into account in the analyses 
Bl to determine the exponents in the main term of Sjy(t). 
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The paper is organized as follows. The asymptotic evaluation of Sjv(£) on fractal lattices is presented in Sec. Q. 
The mathematical techniques involved in the calculation are very similar to those corresponding to the Euclidean case 
and we will only outline the main steps. Details may be found in Ref. In Sec. we compare the asymptotic 
expansion of S'jv(t) with simulation results obtained on the Sierpinski gasket. We end in Sec. |V| with some remarks 
on the quality of the asymptotic approximation. 



II. TERRITORY COVERED BY TV RANDOM WALKERS ON A DETERMINISTIC FRACTAL 

SUBSTRATE 

We will define the multiparticle survival probability T^(t, r) as the probability that site r has not been visited by 
time t by any of the N random walkers that start from the origin site r = at time t = 0. From this definition, the 
following relationship between the average number of distinct sites visited, S'jv(t), and the survival probability can be 
established |,|: 

{l-r w (i,r)} , (4) 

r 

where the sum is over all the sites in the fractal lattice. For independent random walkers, we have T^{t, r) = [r t (r)] W , 
where Tt(r) = Ti(t, r) is the one-particle survival probability. Since we are interested in the behavior of Sjv(i) after a 
large number of steps (thus beyond of the very-short-time regime I), we replace T t (r) and Sjv(i) by their continuum 
approximations: 

/>oo 

S N (t)= [l-r?(r))dV(r) (5) 
Jo 

= V (2D) d f/ 2 t d f/ d ™J N (d f ), (6) 

where we have assumed that the fractal volume V(r) (i.e., the number of lattice sites) of the hypersphere with radius 
r is given by Eq. (||) and that, essentially, Tt(r) only depends on the distance r. In Eq. (||), Jjv(d/) is 

j N (d f )= N l o °° mr- 1 ^-^^. (7) 

In order to evaluate this integral for N 3> 1 it suffices to know T t (r) for large £. For Euclidean lattices 

r t (r) « m ee 1 - iff" exp(-cC) U + fj Kr na ^j (8) 

when £ 1 In this paper we shall consider only those fractals, such as the Sierpinski gaskets, for which 

this equation approximately [see discussion below Eq.(^)] holds. Indeed, we can show that Eq. (||) is reasonable 
for those lattices (e?-dimensional Sierpinski fractals, Given-Mandelbrot curve, one-dimensional lattice, hierarchical 
diamond lattice,. . . ) where the renormalization procedure implemented in Refs. |l7| - |l9[ can be set up. The argument 
is as follows. Let r„ , n = 1, . . . , z, with = |r„ |, be the position of the z nearest neighbors of the site at r = 
in the fractal lattice decimated i times (see Fig. Q), and let h(t,r^') be the probability that, in the time interval 
[0,t], a single diffusing particle that starts at r = is absorbed by any of the traps located at its z nearest neighbors 

placed at the sites • For large values of £ ~ / R, i.e., for relatively short times, the event "the random walker 
arrives for the first time at site r^- 1 " and the event "the random walker arrives for the first time at site r„" are 
(almost) independent so that h(t,r"') is (almost) the sum of the probability [given by 1 - r t (r«)] of each of the 
z individual events, i.e., h(t, r™) ~ z [l ~T t (r^)\. Of course, the two events are not fully independent because 

(i) . (i) 

the random walker could first arrive at site after passing by the site r„ . However, this is very unlikely because 
Tm and rl'' are separated by distances of order of rW so that the fraction of random walkers that, after arriving 
at Yn\ travel to r„ in the short time interval [0,t] is of the order exp(— £"), with £ ~ r^/R 3> 1, because the 
propagator or Green's function is of this order for large values of £ [|l^,p0f. Therefore, it is reasonable to assume that 
1 - r 4 (rM) = h(t 7 r^)/z{l + 0(exp(-£ u )]} for large £. But 1 - h{t,r®) has the form of the right-hand side of Eq. 
(g) with u — d w /(d w — 1), at least for the fractals that we are considering (T^JT^[, so that Eq. (0) for r = rW follows. 
In this discussion we have assumed, in order for the renormalization analysis that leads to T t (rW) to work, that the 



3 








FIG. 1. The 11-generation Sierpinski lattice used in the numerical simulations where the smallest triangles represent for 
8-generation lattices. The sites labelled A t , B t , d and D 4 are nearest neighbors of in the i-th times decimated lattice E^Egl . 
In case I all random walkers start from the origin O. In case II the common origin is chosen randomly from among the sites in 
the shaded area. 



traps were placed at those sites (such as A gi B gi . . .in Fig. |l]) that become the nearest neighbors of the starting site 
after several decimations. But, in the calculation of S^it), the function T t (r) is required for every pair of origin and 
destination sites in the lattice. At this point, we shall assume that the survival probability for r ^ rW and r = rW 
are very similar, i.e., we assume that I\(rW) ~ T T (r) — r r (£) when r/r 1 ^*" = /t 1 / dw — (2D) 1 / 2 ^ for large £, 
with r t (£) given by Eq. (|8|). To the best of our knowledge this problem has not been studied and will require a 
specific and detailed simulation analysis that is not the object of the present work. Nevertheless, previous simulations 
on the Sierpinski gasket of other statistical quantities closely related to F t (r), such as the propagator (or Green's 
function) [ pp| , enable us to affirm with confidence that the parameters c and u remain unchanged over the whole 
Sierpinski lattice whereas the subdominant ones, A, /i, hi,. . . , do not. In the following, we use c — 0.981 |lj],|2(| and 
u = d w / (d w — 1) ~ 1.756 for the two-dimensional Sierpinski gasket (values of c for other fractals can be found in Ref. 
fl7|| ) but we must bear in mind that the actual values of A, \x and h n , n = 1, 2, . . . will likely differ from those obtained 
by renormalization techniques (namely, A = 0.61, fi — 1/2, hi — —0.56) as these latter correspond to the special 
placing of the traps and origins. We must also point out that, as shown below, since the parameters h n , n — 1,2,3..., 
only contribute to the second and higher order series terms of S^it) and since the real value of even the first-order 
asymptotic term is uncertain, the values of these parameters will not be considered here. 

The evaluation of JN(df) now proceeds in analogy with the analysis for the Euclidean lattices |7],^| and we shall 
only quote the final result for Sjv(i) in the fractal case: 

S N (t)^S N (t)(l~A) (9) 

with 

hiN^ df/v 



S N (t) = V (2D) d f/ 2 t d ^[^j (10) 

oo n 

A^f3^2 N s m ) ln ™ ln N > t 11 ) 



n— 1 m— 



and where, up to second order (n = 2), 



S « = -co (12) 



(13) 



4 2) = -(0 - 1) (j2 + t) - {chi ^ (14) 

4 2) = -/i 2 + 09 - (15) 
4 2) = -^(/3-l)^ 2 (16) 
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FIG. 2. Basis vectors used in the simulations to identify the lattice sites of the two-dimensional Sierpinski lattice. The 
labels R, L and C correspond to the three different site types according to the relative position of their neighbors. 



and u> = 7 + In A + /line, 7 ~ 0.577215 is the Euler constant, and /? = df/u = dt(l — l/d w ). 

III. NUMERICAL RESULTS FOR THE SIERPINSKI GASKET 

To check the reliability of the analysis presented in the previous section we carried out simulations of the number 
of distinct sites visited by N random walkers on a two-dimensional Sierpinski lattice with g — 11 generations. This 
means that if we take the length of any side of the smallest triangles (the zeroth decimated triangles) as the unit 
length, then the length of the sides of the triangle that inscribes the lattice (the gth decimated triangle) is 2 9 . Two 
different cases are analyzed: (I) random walkers are initially placed upon the center of the base of the main triangle 
which inscribes the lattice (point O in Fig. Q), and (II) the common starting site is randomly selected. Qualitatively 
and quantitatively, the results are different in case I and case II. The structure of the lattice gives rise to oscillations 
superimposed on the general trend of Sjv(i) in case I. This structure is smeared out in case II by the double average 
over experiments and over starting sites, so that Sjv(i) is now a smooth function. 

We shall use the top vertex of the main triangle as the origin of the lattice and the orthogonal basis {e„, e^,} shown 
in Fig. |^, because it is convenient to divide the lattice into horizontal sets of sites with a fixed v coordinate. The 
positions of the random walkers are updated at every time step by a random selection of the destination sites among 
the neighbors of the site occupied by the walker. In order to do that we have to take into account that the sites of 
the Sierpinski lattice may be classified into three types according to the relative positions of their neighbors: type R, 
L and C as shown in Fig. [| Since the Sierpinski lattice has many holes on every scale, it is not efficient to identify 
the lattice sites and the random walkers positions by two coordinates {v, h} because there are many coordinate pairs 
corresponding to no actual lattice site. Much memory can be saved if the sites are numbered from top to bottom and 
from left to right, so that the top vertex of the main triangle is site number 1 and the right vertex is site number 
3(3 ff + l)/2. 

A. Case I 

First, we will discuss the simulation results for the territory covered by N random walkers placed initially at site 
O in Fig. |. In order to compare with the zeroth- and first-order asymptotic expression, cf. Eq. (^), we must know 
the values of Vo, D, c, u, A and /i. In Fig. ||(a) we have plotted the fractal volume of a circle of radius r, V(r), 
centered upon the privileged site O. The observed log-periodic structure is a consequence of the empty and filled 
triangular areas that repeat periodically as r is increased, but the general trend is well represented by a term of the 
form V r d f with V Q ~ 3.0 ± 0.1. In Fig. |(b), in which the quotient V(r)/r d f is plotted versus log 2 r, one clearly 
sees the log-periodic oscillations of Vq. As our theory assumes a constant value for Vo, we take the average value 
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FIG. 3. The fractal volume V(r) of the Sierpinski lattice as a function of r for case I. (a) log 2 V(r) versus log 2 (r); (b) 
Vo=V(r)/rf versus log 2 (r). 



over the last period (from maximum to maximum), Vq — 2.93, as a reasonable criterion for comparison with the 
simulation results. In order to find the diffusion coefficient D of a random walker starting at O, we performed 10 6 
simulations up to t = 400. The linear numerical fit between t = 50 and t = 400 gives d w ~ 2.32 (the exact value is 
d w — In 5/ In 2 ~ 2.322) and 2D ~ 1.05. Numerical fits using other time intervals (excluding short times, of course) 
lead to similar values, and we take 2D ~ 1.05 ± 0.02 as a reliable estimate. For the parameters c and u we take 
the values (see Sec. |j) 0.981 and d w /(d w — 1) = 1.756, respectively. As discussed in Sec. B the values of [i and A 
are much less certain and we will use here two pairs of values: those obtained by renormalization, i.e., fi = 1/2 and 
A = 0.61, and these same values increased by one, i.e., /x = 3/2 and A = 1.61. Of course, this last pair of parameters 
are arbitrary (other values could also be used) and are mainly given to show the relevance of the corrective terms. 

Simulation results for SV(i) until t — 1000 are shown in Fig. [| for N = 10 3 and N — 10 6 . Overall good agreement 
is obtained in the comparison with the theoretical prediction of Eq. (^|), especially when the values fi — 3/2 and 
A = 1.61 are used. This last finding should only be taken as the manifestation of the importance of the corrective 
terms which can lead to such dramatic changes and improvements in S'jv'(t) after modifying some of the subdominant 
factors of the survival probability. Of course, additional independent study will be necessary to check the form of r t (r) 
given by Eq. (|^) and to find out whether the values for fi and A considered here are good estimates of the real values 
pl[ . Nevertheless, it should be noted that the decrease in 1 — T t (r) when averaging over the whole lattice with respect 
to its renormalization value, as is implied by the corresponding increment of jj, (/x = 1/2 — ► [i = 3/2), is analogous 
to the decrease of the propagator when this same averaging is performed. Given that the two statistical quantities 
(propagator and survival probability) are closely related, one is inclined to accept that, at least, the proposed increment 
in the value of fi captures the right tendency. The subsequent improvement in the prediction of £jv(i) supports this 
supposition. 
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FIG. 4. The number of distinct sites visited on the Sierpinski lattice until t = 1000 by N = 10 3 (circles) and N = 10 6 
(squares) random walkers. All random walkers start from the origin O shown in Fig. [j]. The dashed line is the zeroth-order 
approximation, the dot-dashed line corresponds to the first-order approximation with n — 1/2 and A — 0.61 and the dotted 
line is the first-order approximation with /i = 3/2 and A = 1.61. The solid line also corresponds to this last approximation but 
using for Vo the numerical values of the last oscillation shown in Fig. ^|(b) 



The theoretical expression was not able to give a perfect account of the log-periodic oscillations superimposed on 
the general trend of Sjv(i) shown in Fig. [|. Notice the self-affinity of the numerical Sjv(i) plot (the analytical lines 
are obviously self-affine): if the segment of horizontal axis of Fig. [| between t — and t — 200 is expanded by a 
factor 5 and the corresponding segment of the vertical axis is expanded by a factor 3 = 5 H (H = df/d w — d s /2), we 
get a figure indistinguishable from Fig. [|. The origin of these scaling factors is clear: the enlargement of the sides of 
the Sierpinski gasket by a factor 2 implies that its fractal volume increases by a factor 3 and the time that a random 
walker takes to traverse it increases by 5 || . 

The structure of Sjy(t) is more clearly perceived in Fig. || where the quotient between the theoretical prediction 
and the simulation results is plotted. It is remarkable how relatively poor is the performance of the zeroth-order 
approximation (or main asymptotic term) in predicting the value of Sjy(t): it accounts for hardly eighty per cent of 
Sjv(i) for values of N as large as 10 6 . However, the inclusion of the first corrective asymptotic term (especially for 
some suitable selections of the subdominant parameters A and n) leads to a noticeable improvement. The log-periodic 
structure is observed both for N = 10 3 and TV = 10 6 but in the latter case this structure is richer and strikingly 
similar to that of V(r) as shown in Fig. ||(b). We attribute this fact to a better mapping of the lattice structure as 
more and more random walkers are involved in the exploration. We have plotted the solid line in Fig. |^ with the aim 
of showing to what extent the oscillatory behavior of Sjv(i) as shown in Fig. |^ can be interpreted as a consequence 
of the oscillatory behavior of Vo shown in Fig. |](b). The line is generated in the same way as the dotted line, i.e., 
by means of the first-order approximation of Eq. (||) with A — 1.61 and /x = 3/2, but, instead of using the averaged 
value Vo = 2.93 (as in the dotted line), we use the actual oscillatory value of Vq taken from the last oscillation (from 
maximum to maximum) shown in Fig. ||(b). The way in which the solid line runs alongside the simulation results 
supports the view that the log-oscillatory behavior of SW(£) mainly comes from the log-oscillatory behavior of Vq. 

B. Case II 

We also study the effect on Sjy(t) of choosing other lattice sites, besides the point O, as starting sites for the random 
walkers. To this end we performed simulations where all the N random walkers start on a site randomly selected 
within the shaded area of Fig. [j] in order to avoid the finite size effects. As expected, the fractal volume V(r) and 
the average number of distinct sites visited Sjsf(t) are smooth functions in this case. An estimate of Vq by numerical 
fitting gives Vq ~ 3.6. The analysis of the simulation results (10 4 runs for 10 3 randomly selected starting sites) for the 
mean-square displacement of a single random walker is compatible with 2D « 0.8 when the fit is carried out inside 
the time interval (t — 50, t — 400). Simulation results for Siy(t) (five runs for 10 3 randomly selected starting sites) 
until t = 200 for N — 1024 are shown in Fig. |j. They are compared with the theoretical prediction of the zeroth- 
and first-order approximations of Eq. (^J) for fi = 3/2 and A = 1.61, and with the corresponding simulation results 
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FIG. 5. The ratio 1Z between theoretical and numerical values of S'jv(i) for the two-dimensional Sierpinski lattice with 
N — 10 3 (dashed line) and N = 10 6 (solid line). The asymptotic approximations considered are (from bottom to top) the 
zeroth-order approximation, the first-order approximation with fj, — 1/2 and A — 0.61, and the first-order approximation with 
fj, = 3/2 and A = 1.61. 



when the origin is at O (case II) . Again, we find a relatively poor performance of the zeroth-order approximation, as 
well as substantial improvement when the fisrt-order approximation is used, although there is still room for further 
enhancement. It should be noted that the performance of the two asymptotic approximations is completely analogous 
to that obtained for Euclidean lattices p]||- In these Euclidean media we found that the second-order asymptotic 
approximation gives rise to a significant improvement in the estimate of <S/v(£) even for relatively small values of N. 
It thus seems natural to conjecture that the same will occur for the Sierpinski gasket, although definitive confirmation 
of this guess must wait until reliable values for A, /i and h\ are calculated. 

Finally, in Fig. |?] we show the dependence on N of Sjv(i) for case II. We have plotted two first-order asymptotic 
curves: for the first curve we take the usual values fj, = 3/2 and A = 1.61, and the new values [i — 1.75 and A = 1.75 
are used for the second curve. Again, one sees the great importance of the asymptotic corrective terms as they 
substantially improve the zeroth-order (main term) asymptotic prediction. We have used the new pair of parameters 
simply as another example to illustrate the gross effect of the subdominant factors of the survival probability T t (r) 
on the theoretical prediction of Sj^(t). The excellent agreement reached with /j, = 1.75 and A = 1.75 should not, 
however, be considered as an indication that they are the correct subdominant parameters of the survival probability 

©• 

Notice that the simulation results and theoretical predictions for SN(t) differ very little from case I to case II despite 
Vb and D being clearly different in the two situations. The reason for this coincidence is that the quantity Vo(2D) df / 2 
that appears in the amplitude of the main asymptotic term of S]y(t) is almost invariant with respect to translations 
of the origin site: its value was approximately 3.05 when the origin is placed at O, and 3.02 when the origin sites were 
randomly selected. 

IV. SUMMMARY 

We addressed the problem of calculating the average number 5at(^) of distinct sites visited (or territory covered) 
by N independent random walkers that diffuse on deterministic fractal lattices. The validity of the main result of 
this paper, Eq. (^|), is based upon two conditions: (a) the volume (number of sites) of the fractal substrate grows as 
V{r) = Vor d f and (b) the asymptotic form of the survival probability, T t (r), for large r and small t has essentially 
the same form as that corresponding to the Euclidean lattice case. The mathematical method used to derive such a 
result had already been successfully applied to Euclidean lattices 0|| and the fractal case is a fairly straightforward 
generalization when the previous conditions are fulfilled. In order to check the goodness of the approximation, we 
carried out numerical simulations on a standard deterministic substrate (the two-dimensional Sierpinski gasket) ob- 
taining reasonable agreement with the theoretical results, especially when theoretical first-order asymptotic corrective 
terms are considered. The performance of the theoretical expressions discussed closely resembles that attained for 
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FIG. 6. The average number of distinct sites visited in the Sierpinski lattice for N = 2 8 (filled circles) and N = 2 13 (filled 
squares) until t = 200 when the common origin is randomly selected from the shaded area of Fig. [j]. The open symbols represent 
the corresponding values when the origin is the point O. The solid lines are the zeroth-order approximation and the dashed 
lines are the first-order approximation with /i = 3/2 and A = 1.61 for, from top to bottom, N = 2 13 and N = 2 s . 




FIG. 7. Dependence on N of the fractal territory Sjv(t) explored by N random walkers by time t = 100 (circles) and t = 200 
(squares) when the common origin is randomly selected from the shaded area of Fig. |l|. The solid lines are the zeroth-order 
approximation, and the dashed lines [dot-dashed lines] correspond to the first-order approximation with /i = 3/2 and A — 1.61 
[fj, = 1.75 and A = 1.75] for, from top to bottom, t = 200 and t = 100. 
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Euclidean lattices. However, a more definitive check of the theoretical expressions for <SW(£) that include corrective 
terms is hindered by the uncertainty in the value of the subdominant parameters A, fi, hi, . . .that appear in the 
survival probability T t (r). Unfortunately, this function cannot be completely determined by the rigorous renormal- 
ization scheme mentioned in Sec. [n], so that the definitive determination of its subdominant terms by numerical (or 
other analytical procedures) is a problem for future work which will surely be beset by with the technical difficulties 
associated with the identification of these faint subdominant terms [EOfl . 
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